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REMARKS 

The amendment filed August 22, 2005 was not entered, so the present changes are based 
on the claims as they appear in the Amendment filed November 26, 2004. Claims 51-72 remain 
active in the case. Claim 73 has been cancelled as redundant in view of amended claim 5 1 . 

Support for classifying a test substance in one or more functional categories based on 
reference to a data library formed from known compounds is found on page 30, lines 11-12; 
page 31, lines 7-15; page 32, Table 1 and page 33, lines 1-17. 

Serum-free culture medium is supported on page 26, lines 4-6 ("The neurons are 
grown. . .in serum-free media. The cells are cultured in defined media. . ."). 

An intervening layer acting as a high-impedance seal is supported in the specification at 
page 11, lines 6 - 12. 

Measurement of changes in one or more characteristics selected from after potential, time 
to cessation of activity, frequency, amplitude, shape, spike rate, or time constant_are described on 
page 40, lines 14-31 and this subject matter was also in claim 53. 

Deconvoluting changes to identify one or more ion channels affected by the test 
substance is supported on page 7, lines 13-28 and page 10, lines 4-16. Please note we 
specifically define deconvolution as not including spectral analysis on page 13, lines 1-5 . 
Therefore, the term "deconvolution" does not read on Fourier Transform analysis. 

Our initial description of "functional categories" came from the literature (see ref Riley et 
at in the application as well as page 6 lines 17-18) and it is a method of classifying related genes 
and their gene products and thus the pathways they enable. We have shown differences in 
Action Potential peak shape with the data included in the application, see Figure 5. 

We have since demonstrated the application in our paper "Toxin Detection based on 
Action Potential Peak Shape Analysis using a Realistic Mathematical Model of Differentiated 
NG106-15 Cells" that has been accepted and is in press in the journal Biosensors and 
Bioelectronics (Mohan, et al., submitted herewith). 

Table 1 on page 32 of the application also further illustrates our ability to group the 
effects of different compounds into different "functional categories". 
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As discussed on page 40 lines 21-32 and page 41 lines 1-4 of the application, 
deconvolution of an action potential implies that some prior knowledge of the factors 
contributing to the shapes are known. It is well known that an action potential requires the 
contributions of ion fluxes through at least three distinct ion channels. Varying the contributions 
of these ion fluxes from these channels has been shown to very the shape of the action potential, 
and this is illustrated in Figures 5 and 6A. of the application. Thus, to deconvolute an action 
potential into its components and relate them to pathways or functional categories inside a cell 
requires some prior knowledge of the biology of the system (i.e., a data base of biological test 
data on the subject cell culture) in order to perform the deconvolution step and identify cell 
pathways. The Mohan, et al. paper illustrates some of our efforts in this area and an 
implementation of the algorithm shown in Figure 10 of the patent application. This algorithm 
has been used to deconvolute one of the examples given in Figure 5 of the application and is 
shown in Figure A, below. 

Control Ricin 1 mM 




Time (ms) Time (ms) 

Figure A. Result of parameter fitting to action potential shapes before (left panel) and after (right 
panel) the application of 1 mM Ricin. Action potentials were recorded from differentiated 
NG108-15 cells using whole-cell patch clamp method. Action potentials was fitted (dotted lines) 
with a computer model of NG108-15 cells established based on the measurement of ion channel 
parameters. 
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Table 1. Effect of 1 mM Ricin on the action potential parameters. Action potential parameters 
were obtained by fitting the parameters of our NG108-15 cell model to the experimental data. 1 
mM Ricin caused a significant slow-down of the sodium channel kinetics (parameter A) and a 
decrease of sodium conductances in the membrane (parameter g). The effect of Ricin on the 
other ion channels were not so profound, partly because potassium and calcium channels are 
activated at higher membrane potentials, so their contribution to the Ricin-modified action 
potential is less (the amplitude of the action potential is smaller, barely reaching 0 mV). 



This type of analysis is very different than a spectral analysis "i.e. a Fortier Transform" 
which arbitrarily assigns combinations of sign waves of different frequencies and amplitudes, to 
represent the area under a curve. The assignment of the different sign waves is arbitrary and, 
although it gives a unique signature to a signal, there is no way to relate it to biological 
information, except through an empirical approach to the testing of drugs or other compounds as 
well as toxic chemicals. 

Claims 51-67 and 70 - 73 are rejected under 35 USC 103(a) over Borkholder et al. in 
view of Georger, Jr. et al. 

Borkholder et al. does not specifically disclose an intervening layer, as claimed. 
Furthermore, it does not disclose or suggest serum-free media. Finally, it does not inherently 
disclose a deconvolution analysis based on a data library, as claimed. Instead it mentions that a 
"pattern of spectral changes provides a unique signature for the compound. Matching the 
spectral change pattern of a test compound with the spectral change pattern for reference 
compounds provides a useful method for characterization of channel modulating agents." (Col. 
4, lines 15-20, emphasis added). This passage describes using the observed changes as a 
fingerprinting method to associate certain channel-modulating compounds with the observed 
changes, but that is not the same kind of process as deconvoluting a signal to identify the 
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contributions of different ion channels to the signal and then to assigne the test compound to one 
or more a functional categories. 

Georger Jr. et al. teaches a method of making biosensors using photolithography 
placement techniques. It never explicitly states that a high-impedance seal is created or that it 
would enable deconvolution. For instance, col. 7, lines 58-63 does not relate to deconvolution, 
but at most implies a high-impedance seal. Col. 10, lines 58-62 describes manual placement, 
whereas we deposit the cells by cell culture. At col. 13, lines 19-24, the cell culture patterning 
method is to use fewer steps. The discussions at col 16, line 45 and at col. 10, lines 50 - 65 relate 
to positioning. 

Finally, Georger Jr. et al. uses serum (col. 17, lines 60 - 65) and states at col. 18, lines 9 - 
10 that "adhesion is unaffected by serum." Especially where high-throughput screening is 
involved, the absence of serum is vital. 

Claims 51-73 are rejected under 35 USC 12, first paragraph, for written description 
problems. Support for the high-impedance intervening layer allowing deconvolution as claimed 
is found on page 11, lines 6 - 12. 

Claim 51 is rejected under 35 USC 12, first paragraph, on enablement grounds. A 
positive deconvolution step is now recited. 

Claims 51 - 73 are rejected under 35 USC 12, second paragraph, as being indefinite for 
failure to positively recite a deconvolution step. This has been address by amendment. 

Claim 57 has been clarified with respect to the position of the insulator "surrounding" the 
electrode. 

The dependency of claim 61 has been corrected. 

Claim 67 does add a limitation over claim 66 because "known function" and "unknown 
function" are not the only two possibilities for DNA; some DNA has no function at all. 
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Abstract 

//■ \\ . . 

The NG 108- 15 neuroblastoma/glioma hybrid cell line has been frequently usedTor toxin detection, pharmaceutical screening and as a whole-cell 
biosensor. However, detailed analysis of its action potentials during toxin or drug administration has not been accomplished previously using patch 
clamp electrophysiology. In order to explore the possibility of identifying toxins based on their effect on the shape of intracellularly or extracellularly 
detected action potentials, we created a computer model of the action^potential generation of this cell type. To generate the experimental data to 
validate the model, voltage dependent sodium, potassium and high-threshold.calcium currents, as well as action potentials, were recorded from 
NG 108-15 cells with conventional whole-cell patch-clamp methods. Based on the classic Hodgkin-Huxley formalism and the linear thermodynamic 
description of the rate constants, ion-channel parameters were estirnated;using an automatic fitting method. Utilizing the established parameters, 
action potentials were generated in the model and were optimized to represent the actual recorded action potentials to establish baseline conditions. 
To demonstrate the applicability of the method for toxin detection^anddiscrimination, the effect of tetrodotoxin (a sodium channel blocker) and 
tefluthrin (a pyrethroid that is a sodium channel opener) were studied. The two toxins affected the shape of the action potentials differently and their 
respective effects were identified based on the changes in the fitted parameters. Our results represent one of the first steps to establish a complex 
model of NG 108- 15 cells for quantitative toxin detection based oh action potential shape analysis of the experimental results. 
© 2005 Elsevier B.V. All rights reserved. // ^ 

Keywords: Action potential shape analysis; Toxin detection; tyG 108-1 5; Computer simulation; Linear thermodynamic model; Hodgkin-Huxley model 



1. Introduction / /? n 

:■ yf :/ 

In the areas of environmental protection; toxicology and drug 
development there are increasing demands for high- throughput 
functional screening methods (Rogers,^ 1995; .Paddle, 1996; 
Ohlstein et al., 2000; Heck et al., 200lf Grosjton, 2002; Tzoris 
et al., 2002). For monitoring of the environment, whole-cell 
biosensors could be more effective than physico-chemical meth- 
ods to assess the global toxicity of the wide variety of chemi- 
cals that are possible pollutants (Evans/et al., 1986; Rogers, 
1995; Bousse, 1996; Paddle, 1996; Naessens and Tran-Minh, 
1998a,b). Whole-cell biosensors are also able to give functional 
information about the effect of chemicals, have the ability to de- 
tect unknown compounds and continuous monitoring of external 
conditions is possible as theyxan be made small enough to allow 

* Corresponding author. Tel:: +1 864 710 8472; fax; +1 407 882 1 156. 
E-mail address: jhickman@mail.ucf.edu (J.J. Hickman). 



field applications (Bousse, 1996). Another benefit of whole cells 
in environmental applications is that they allow the measurement 
of the total bioavailability of a given pollutant rather than its free 
form (Bousse, 1996; Naessens and Tran-Minh, 1998a,b; Phiip 
et al., 2003). Similarly, in safety pharmacology, the side effect 
spectrum of a given compound is not known, thus, the appli- 
cation of complex functional tests at the whole-organ ism- level 
are necessary (Jorkasky, 1998; Kinter and Valentin, 2002) and 
could be addressed with this technique. Moreover, the availabil- 
ity of genomic information significantly increased the number 
of potential targets available for drug discovery and new meth- 
ods are necessary for high- throughput funptional screening for 
target validation (Ohlstein et al., 2000; Croston, 2002). 

Recently, the application of whole-cell biosensors for 
toxin detection and drug screening has become more readily 
accepted (Bousse, 1 996; Bentleya et al., 200 1 ; Baeumner, 2003) 
as it has many benefits compared to traditional methods of 
evaluation. Several techniques have been developed to quantify 
the physiological changes induced by chemicals in whole-cell 
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sa biosensors (Bousse, 1996; Bentleya et al., 2001). One of 

59 these techniques, which is frequently used for monitoring the 

eo physiological state/activity of excitable cells, is multi-electrode 

61 extracellular recording of membrane potential (Bousse, 1996; 

Gross et al., 1997; Denyer et al, 1998; Jung et al, 1998; 

es Offenhausser and Knoll, 2001; Krause et al., 2000; Stett et 

64 al., 2003). The high-throughput or long-term application of 

es extracellular recording is much preferred over intracellular 

ee action potential recording for many applications because the 

67 use of an intra-cellular or patch clamp electrode limits the life 

es of the cell to a few hours as does the use of voltage sensitive 

69 dyes (all dyes reported to date, to a greater or lesser extent, are 

70 toxic to cells) (Mason, 1993; Chiappalone et al., 2003). 

71 Action potential generation and the shape of the action po- 

72 tential depends on the status of several ion channels located in 

73 a cell's membrane, which are regulated by receptors and intra- 

74 cellular messenger systems (Gross et al., 1995, 1997; Morefield 

75 et al., 2000). Changes in the extracellular (receptor activation) 

76 or intracellular environment (gene expression), in many cases, 

77 can be reflected in an alteration of spontaneous firing proper- 

78 ties such as the frequency and firing pattern (Gross et al., 1997; 

7 9 Amigo et al., 2003; Chiappalone et al., 2003; Xia et al., 2003) 
so of excitable cells and also in changes in the shape of their action 
si potentials (Clark et al, 1993; Muraki et al., 1994; Akay et al., 
82 1998; Nygren et al., 1998; Djouhri and Lawson, 1999). There 
es are many examples indicating that the shape of action potent 
&4 tial depends on the extracellular and intracellular environment 

85 of the cells. Sodium (Spencer et al., 2001), potassium (Claris 

86 et al., 1993; Marti n-Caraballo and Greer, 2000) and calcium 
channel modulators (Ahmed et al., 1993; van Soest and^Khs,, 

as 1998), as well as several toxins and various pathological condi^ 

89 tions (Muraki et al., 1994; Shaw and Rudy, 1997; Akay efaL, 

90 1998; Djouhri and Lawson, 1999) have already been' shown to 

91 affect the shape of action potentials. However, action potential 

92 shape analysis for high-throughput screening applicVtions^such 

93 as toxin detection or drug screening, has not yet t>een "developed. 

94 Two primary reasons this type of analysis has been lacking to 

95 date is that it is difficult to obtain high fidelity-recordings from 

96 chip-based extracellular electrodes and from^hehackx)f models 

97 to adequately analyze the signals. "'"^C.^ 

98 Several mathematical models have been developed to de- 

99 scribe the electrical properties and the -process of action po- 

100 tential generation in excitable cells (Agin,"1972; Cohen, 1976; 

101 Otten and Scheepstra, 1995; Dokos and Loyell, 1996; Weiss, 

102 1996; Shevtsova et al., 2003). The' 1 most ^widely used is the 

103 Hodgkin-Huxley formalism where ion channel activation and 

104 inactivation is described using voltage dependent activation and 

105 inactivation gates (Weiss, 1996). In theofiginal model the volt- 

106 age and time dependence of the gates was given utilizing rate 

107 constants, which were taken as an empirical function of the mem- 

108 brane potential (Hodgkin and Huxley, 1 952). However, in lieu 

109 of using empirical functions, it is also possible to deduce the 
no functional form of the voltage dependence of the rate constants 
m from thermodynamics (Weiss, 1 996; Destexhe and Huguenard, 

112 2000). The applications of these models to extracellular record- 

113 ings from a suitable excitable cell population, in combination 

114 with the proper models, would be a logical next step in adapt- 



ing this technology to high throughput toxin detection and drug n 5 

discovery. 116 

The NG 1 08- 1 5 hybrid cell line, which was created by merg- m 

ing mouse neuroblastoma and rat glioma cells, has been widely us 

used in in vitro experiments as a substitute for primary-cultured ns 

neurons (Hu et al., 1 997; Doebler, 2000; Tojima et al, 2000). The 120 

neuronal functions and features of differentiated NG 1 08- 1 5 cells 121 

have been well characterized, e.g. the presence of a wide range of 122 

voltage dependent and transmitter activated membrane currents 123 

have been detected as well as second messengers and enzymes 124 

normally found in primary neurons (Schmitt and Meves, 1 995; 125 

Lukyanetz, 19^8; Ma et al, 1998; Tojima et al., 2000). NG108- 126 

1 5 cells are widely used in pharmacology (Hu et al., 1997) and 127 

also as a whole-cell biosensor for toxin detection (Ma et al., 12s 

1998). Oneof the ) distinctive features of the NG108-15 cell line, 129 

which makes Jtiideal for whole-cell biosensor applications, is 130 

that the^cells do not form synaptic connections, thus network 131 

activity does riot influence single cell data (Ma et al., 1 999). 132 

In this study we created a computer model of the action po- 133 

tentialgeneration of an NG1 08- 1 5 cell based on voltage-clamp «4 

and^cuffent clamp electrophysiological recordings. Using this 135 

model we developed and demonstrated the applicability of ac- 136 

^tic^n potential shape analysis as a method for toxin detection and 137 

monitoring the physiological state of excitable cells. 133 

2i Methods 139 

^ 2.7. Surface chemistry mo 

NG108-15 cells were plated on N-l[3-(trimethoxysilyl) «i 

propyl] diethylenetriamine (DETA) coated glass coverslips 142 

(22 mm x 22 mm, Thomas Scientific). The DETA coated cover- 143 

slips were prepared according to published protocols (Schaffher 144 

et al., 1995). In brief, glass coverslips were cleaned using us 

HCl/methanol (1:1) followed by a concentrated H2SO4 soak we 

for 30 min followed by a water rinse. The coverslips were then 147 
boiled in deionized water followed by a rinse with acetone and 

then oven dried. The DETA films were formed by the reac- 14 9 

tion of the cleaned surfaces with a 0.1% (v/v) mixture of the 150 

organosilane in toluene. The DETA cover glasses were heated 151 

to just below the boiling point of toluene, rinsed with toluene; re- 152 

heated to just below the boiling temperature again and then oven 153 

dried. 154 



2.2. Culture ofNG108-15 cells 

The NG108-15 cell line (passage number 16) was obtained 
from Dr. M.W. Nirenberg (NIH). The NG108-15 cells were cul- 
tured according to published protocols (Higashida et al., 1986; 
Ma et al., 1998). Briefly, the cell stock was grown in T-25 and T- 
75 flasks in 90% Dulbecco's modified Eagle's medium (DMEM, 
GIBCO) supplemented with 10% fetal bovine serum and HAT 
supplement (GIBCO, lOOx) at 37 °C with 10% C0 2 . Differen- 
tiation was induced by plating the cells in a serum-free defined 
medium (DMEM + N2 supplement, GIBCO) in 35 mm culture 
dishes at a density of 40,000 cells/dish. 
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165 2.3. Electrophysiological recordings 

166 Whole-ceil patch clamp recordings were performed in a 
16? recording chamber on the stage of a Zeiss Axioscope 2 FS 

168 Plus upright microscope. The chamber was continuously per- 

169 fused (2 ml/min) with the extracellular solution. The compo- 

170 sition of the extracellular solution for the recording of action 

171 potentials was (in mM): NaCI 140, KC1 3.5, MgCl 2 2, CaCI 2 

172 2, glucose 10, HEPES 10. For the recording of potassium cur- 

173 rents 1 fiM tetrodotoxin (TTX) was added to the extracellular 

174 solution. To minimize space-clamp errors, sodium currents were 
us recorded in a 'decreased sodium' extracellular solution contain- 

176 ing (in mM): NaCI 50, TEA-C1 100, CsCl 5, CaCl 2 1, CoCl 2 

177 1, MgCl 2 1, glucose 10, HEPES 10. For the recording of cal- 

178 cium currents, sodium and potassium channels were blocked 

179 with Cs, TEA and TTX. The extracellular solution composi- 

180 tion for the measurement of calcium currents was (in mM): 

181 NaCI 1 00, TEA-C1 30, CaCl 2 1 0, MgCl 2 2, glucose 1 0, HEPES 

182 10, TTX 0.001. The pH was adjusted to 7.3 and the osmolar- 

183 ity was 320 mOsm. The intracellular solutions composition for 

184 recording the action potentials and for potassium channel mea- I 

185 surements was (in mM): kgluconate 130, MgCl 2 2, EGTA 1, 

186 HEPES 15, ATP 5, for sodium channels was CsF 130, NaCI// 

187 10, TEA-C1 10, MgCl 2 2, EGTA 1, HEPES 10, ATP 5 and^;: 

188 for calcium channels was: CsCl 120, TEA-C1 20, MgCl 2 2, 

189 EGTA 1 , HEPES 1 0, ATP 5 (pH 7.2; osmolarity = 280 mOsm);C\ 

190 For selecting L-type calcium channels, 1 yM a>CT x GVI^was ) 

191 used. /) " / 

192 Patch pipettes (4-6 M £2 resistance) were prepared frorn^ 

193 borosilicate glass (BF1 50-86-10; Sutter, Novato, C/(j-wifa£ 

194 Sutter P97 pipette puller. Voltage clamp and current clamp exper- 
ts iments were performed with a Multiclamp 7 00 A (Xxoh>Union 

196 City, CA) amplifier. Signals were filtered at 2 kHz-and digitized 

197 at 20 kHz with an Axon Digidata 1 322 A interface. Data record- 

198 ing and analysis was performed using pClamp 8 "(Axon) soft- 

199 ware. Sodium and potassium currents were measured in voltage 

200 clamp mode using 1 0 mV voltage steps from a - 85 mV holding 

201 potential. To record high-threshold calcium currenfs/ a -40 mV 

202 holding was used. Whole cell capacitance 'and series resistance 

203 was compensated and a p/6 protocol was i us^d/The access re- 

204 sistance was less than 22 M£2. Action potentials were measured 

205 in current-clamp mode using 1 s depolarizing current injections. 
Me Data was saved in text-format and imported into MATLAB for 

207 further analysis. 

208 2.4. Simulation of ionic conductances and action potential 

209 generation in NG108-15 cells 



210 
211 



dynamic changes in the membrane-potential were calculated ac- 21a 

cording to: * 219 

^external Aonic 

dt Cm 

The dynamics of the state variables was given as 221 

dmldt = {m oc - m)/r m . Where g Na , g K , £ C aL, ^ Na , V K , V C aL are 222 

constants (maximum conductances of the channels and reversal 223 

potentials, respectively); m, n, h, e are the state variables, moo, 224 

n oo y hoo , ^00 are the steady- state values of the state variables and 225 

the r-s are their voltage-dependent time-constants. The voltage- 226 

dependence/bfvthe steady-state state parameters and the time 227 
constants were/given using the general thermodynamic formal- 
ism (using :the state- variable m as an example): 



1 



moo = 



and 



.i^tf>-zF/RT{V m -v l/2 ) 



JJ 



^exp(^/ RT )^( v m-v 1/2 ) C0S ^( z/ r/2RT(V^ - V 1/2 )) 



The classic Hodgkin-Huxley model (Hodgkin and Huxley, 
1952) was used for the description of the ion channel currents, 

212 but instead of the original empirical description of the rate con- 

213 stant, the thermodynamic approach (Weiss et al., 1 995 ; Destexhe 

214 and Huguenard, 2000);was applied. Briefly, the total ionic mem- 

215 brane current was described as: 

216 / ionic = / Na + / K + / Ca + /! = £ Na m 3 /z(V - V Na ) 

+gKn\V - V K ) + g C zte\ V - V CaL ) + g,(V - Vi) 



228 
229 



where.z, V\n,A and £ are fitting parameters and V m represents 233 

the membrane potential. 234 

\>, '^As it can be seen from these equations V1/2 corresponds to 235 

the half activation/inactivation potential of the channel and A is 236 

linearly related to the activation or inactivation time-constant. 237 

]The meanings of z and £ are not as obvious: z is related to 233 

the number of moving charges during the opening or closing 239 

of the channel; whereas £ describes the asymmetric position of 240 

the moving charge in the cell membrane. Sodium, potassium 241 

and calcium channel mediated currents, which were recorded in 242 

voltage-clamp mode at different membrane potentials (10 mV 243 

increments, —40 mV, +30 mV range), were fitted in one step us- 244 

ing the above described model, with the corresponding ion chan- 245 

nels included, using the built-in routines (fminunc) of MATLAB 246 

through a custom-made graphical interface. Parameters obtained 247 

from different cells were averaged (/z = 4-6) and considered as 24s 

initial values for the action potential modeling. Simulated ac- 249 

tion potentials were fitted to the experimental data using built-in 250 

functions in MATLAB (fminunc). Fminunc is used to find a 251 

minimum of a scalar function (the error function, see Section 3) 252 

of several variables, starting at an initial estimate. This method 253 

is generally referred to as unconstrained nonlinear optimization. 254 

Because it is finding only local minimums, it is very important 255 

to start the optimization as close to the final result as it is possi- 256 

ble. In our case, the averaged ion channel parameters obtained 257 

from voltage clamp experiments served the initial values for the 25a 

parameter estimations. 259 

3. Results 260 

The NG108-15 cells completed their differentiation process 261 

and a neuronal phenotype was obtained by day 1 0 in vitro (DIV) 262 

in the defined, serum-free medium (Fig. 1A). Electrophysio- 263 

logical experiments were performed on the differentiated cells 264 

between day 10 and 14 in vitro. All of the cells investigated 265 

showed pronounced sodium, potassium and calcium currents in 266 
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4 6 8 
Time (ms) 

Fig. 1. Estimation of ion channel parameters from the voltage-clamp experiments. (A) Phase-contrast picture of an NG 108-1 5 cell with a patch-clamp electrode 
attached to the cell (40 x objective, scale bar =25 jim). (B) Sodium currents recorded at different membrane potentials of — 10, 0, 10, 20 mV in voltage-clamp mode 
(solid line) and the results of the parameter fitting using the Hodgkin-Huxley model and the linear thermodynamic formalism (dotted line). (C) Potassium currents 
recorded at membrane potentials of 0, 10, 20 and 30 mV (solid line) and the fitted curves using the^rhoder(dotted line). 



267 the voltage-clamp experiments. Most of the cells fired one sin- 

268 gle action potential upon depolarization in current-clamp mode, 

269 whereas about 1 0% of the cells were able to fire multiple action 

270 potentials. Only a very small minority of the cells (about 5%) 

271 were spontaneously active. 

272 3.1. Extracting ion-channel parameters from the 

2?z voltage-clamp experiments using the linear thermodynamic 

274 description 

275 Signals from sodium, potassium and high-threshold (L-type) 

276 calcium channels were recorded in voltage-clamp mode usihg;;^ 

277 the Axon's pClamp 8 program with standard protocols (Fig/lB 

278 and C). The data were saved in ASCII format and imported , 

279 into the MATLAB program. A graphical interface was created^ 

280 to fit the mathematical model to the experimental data and 40 

281 visualize the results. To quantify the difference between the fitted 

282 curves and the recorded data the following error- functions%ere 

283 implemented: X // 



3.2. Action-potential shape analysis 



Maximum error: £ Ma x ~ Max(Abs(/?(f„) ^S^Wf where 
R(t n ) is the recorded value and S(t n ) is the simulated data 



at time t n 



if 



Least square: E Lsqua re = - S('„)) 2 . 

Weighted least square: £wLsquare = ££square f rt <30ms 

and £wLsquare = 5 X ^square if t R ^'30 msV 



After several trials it was concluded mat simulations using the 
weighted least square error function gave the most satisfactory 
results because the other error functions occasionally obtained 
a non-inactivating sodium-current component in the simulated 
data. Curves were fitted after anMnitialyO. 1 ms delay to elimi- 
nate the effect of experimental artifacts. In some simulations the 
reversal potential for the ionic conductances was kept constant. 

In general, an excellent fit td^the potassium channel data 
(Fig. IB and C) and an acceptable fit to the sodium and cal- 
cium channel data were achieved. The automatic fitting algo- 
rithm converged in less than ^min running on a Pentium III 
1 GHz personal computer. After averaging the results of 3-10 
experiments the initial parameter values for modeling the action 
potentials were obtained (Table 1). 



Action potentials were evoked with short (2 ms) current in- sw 

jections in current clamp mode either at resting membrane po- 30s 

teritial orata —85 mV holding potential. The following parame- sos 

ters;;were obtained from the patch-clamp recordings and used in so? 

/the modeling: membrane resistance, resting membrane poten- 3os 

/ /tial/ membrane capacitance and injected current. The maximum 309 

conductance of the leakage current (g\) was calculated from the 310 

_ ionic conductances and from the resting membrane potential. 311 

""We used the earlier established, averaged ion-channel parame- 312 

ters as the initial parameters for the action potential fitting. Using 313 

"- voltage dependent sodium, potassium and L-type calcium on- 314 

[} ductances, an excellent fit to the rising and to the initial falling 315 

phase of the action potentials in the NG108-15 cells was ob- 316 

tained (Fig. 2). 317 

We also obtained an excellent fit to the experimental data 31 s 

in the case of the toxin-modified action potentials by modify- 319 

ing only the corresponding sodium-channel parameters (Fig. 2, 320 

Table 2). 321 



4. Discussion 322 

We have developed a mathematical model of the action po- 323 

tential generation in NG108-15 cells and extracted the param- 324 

eters for this model from whole-cell patch clamp experiments. 325 

Utilizing only voltage- sensitive sodium, potassium and L-type 326 

calcium channels we were able to obtain an excellent fit to the 327 

rising as well as to the initial falling phase of the action poten- 32a 

tial . The weak fit to the later falling phase of the action potential 329 

could be due to active conductances, which were not taken into 330 

account in this model. For example, at lease three other cal- 331 

cium channels and also a calcium activated potassium channel 332 

have already been described in NG 108- 15 cells. In this study, we 333 

kept the number of the ion channels modeled to a minimum, due 334 

to the high computational requirements of the parameter fitting 335 

program. 335 

A linear thermodynamic formalism was used to describe the 337 

voltage and time dependence of the ionic conductances, which 33s 

eliminated the need for 'guessing' the function for the voltage- 339 

dependence of the rate constants and the same form (with dif- 340 
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Fig. 2. Effect of toxins on the action potentials of NG108-15 cells. (Upper panel) effect of 0.5 jjlM tetrodotoxin. (Lower panel) effect of 0.5 p,M tefluthrin. (Solid 
line) data recorded in current clamp experiments. (Dotted line) results of the simulation using the mathematical model of the NG108-15 cells. 
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341 ferent parameters) could be used for the characterizatipn of all 

342 the ion channels. {( V\ 

343 One of the limitations of the model is the, high number of 

344 parameters required to describe the ionic conductances. In this 

345 study four parameters were used to characterize each activa- 

346 tion/inactivation gate, thus, for the description ofCthe sodium 

347 channel, a total of 10 parameters was required/ With this high 

348 number of parameters a more detailed study; is needed to prove 

349 the uniqueness and stability of the solutions. '/ 

350 In order to validate the action potential shape analysis for 

351 utilization as a toxin detection method, the effect of two toxins, 

352 tetrodotoxin, a specific sodium channel blocker, and tefluthrin, 

353 a sodium channel opener pyrethroid were analyzed. Changes 

354 in the action potential shape, and;aIso in the fitted ion-channel 

355 parameters caused by the two toxins, \yere measured. An excel- 

356 lent fit to the toxin-modified action, potential shapes by modi- 

357 fying only the appropriate sodium channel parameters has been 

358 achieved. TTX, as expected from a channel blocker, significantly 

359 decreased the maximum sodium conductance (g), but did not af- 

360 feet the voltage dependence of the channel (V\n) (Table 2). Un- 

361 expectedly, TTX affected the activation kinetics of the sodium 

362 channels as well. TTX moderately increased the activation A pa- 

363 rameter by causing a slowing down of the activation of the chan- 

364 nel. One possible explanation could be the existence of different 

365 subpopulations of sodium channels on the NG 108-15 cells with 



different activation and inactivation time constants and different 
TTX sensitivities. 

The major effect of tefluthrin (a channel opener) was to 
slow down (practically remove) the inactivation of the sodium 
channels (A). Tefluthrin also affected the maximum sodium 
channel conductance (g) and the voltage dependence of the 
activation and inactivation of the channel (V1/2), shifting 
appropriate current-voltage (//V) relationships to the left by 
about 15mV. A Similar effect on the voltage dependence of 
the sodium channels was described by Spencer et al. (2001) for 
fenpropathrin, a tefluthrin-Iike pyrethroid. 

In summary, these experiments indicated that we were able to 
decipher and quantify the effects of toxins on ion channels with- 
out actually measuring ion channel currents in voltage-clamp 
experiments. Instead, changes in the shape of action potentials 
measured by patch clamp electrophysiology, combined with a 
validated computer simulation of the cell, were utilized. 

This method could be useful for toxin detection and for clas- 
sification of unknown toxins in environmental protection sce- 
narios or in the detection of biological and chemical warfare 
agents. It could also be extended to functional screening in drug 
development. With the refinement of the model of the cell not 
only those toxins could be identified, which are directly acting 
on ion channels, but also changes in second messenger levels or 
gene expression could also be detected and classified. Moreover, 
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recent advancements in the study of the cell electrode interface 
and microelectrode-fabrication technology indicate that high- 
fidelity extracellular recording of action potential shapes might 
be possible, which opens a new horizon for the high-throughput 
application of our method using extracellular recordings. 

5. Conclusions 

We have demonstrated that toxin effects on ionic membrane 
currents can be quantified based on the measurement of changes 
in action potential shape and a realistic mathematical model of 
action potential generation in NG108-15 cells. Further studies 
are underway to improve the mathematical model, explore the 
applicability of the method for detection and- identification of 
a wider variety of toxins and to extend the technique to enable 
obtaining high-fidelity action potential data with non-invasive, 
high- throughput extracellular electrodes, in order to make high- 
throughput functional toxin detection and drug screening possi- 
ble. 
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